Drift dynamics in microbial communities and the effective community size.

The structure and diversity of all open microbial communities are shaped by individual births, deaths, speciation and immigration events; the precise timings of these events are unknowable and unpredictable. This randomness is manifest as ecological drift in the population dynamics, the importance of which has been a source of debate for decades. There are theoretical reasons to suppose that drift would be imperceptible in large microbial communities, but this is at odds with circumstantial evidence that effects can be seen even in huge, complex communities. To resolve this dichotomy we need to observe dynamics in simple systems where key parameters, like migration, birth and death rates can be directly measured. We monitored the dynamics in the abundance of two genetically modified strains of Escherichia coli, with tuneable growth characteristics, that were mixed and continually fed into 10 identical chemostats. We demonstrated that the effects of demographic (non-environmental) stochasticity are very apparent in the dynamics. However, they do not conform to the most parsimonious and commonly applied mathematical models, where each stochastic event is independent. For these simple models to reproduce the observed dynamics we need to invoke an 'effective community size', which is smaller than the census community size.

The effect of metabolic stress on genome stability of a synthetic biology chassis Escherichia coli K12 strain.

BACKGROUND: Synthetic organism-based biotechnologies are increasingly being proposed for environmental applications, such as in situ sensing. Typically, the novel function of these organisms is delivered by compiling genetic fragments in the genome of a chassis organism. To behave predictably, these chassis are designed with reduced genomes that minimize biological complexity. However, in these proposed applications it is expected that even when contained within a device, organisms will be exposed to fluctuating, often stressful, conditions and it is not clear whether their genomes will retain stability. RESULTS: Here we employed a chemostat design which enabled us to maintained two strains of E. coli K12 under sustained starvation stress: first the reduced genome synthetic biology chassis MDS42 and then, the control parent strain MG1655. We estimated mutation rates and utilised them as indicators of an increase in genome instability. We show that within 24 h the spontaneous mutation rate had increased similarly in both strains, destabilizing the genomes. High rates were maintained for the duration of the experiment. Growth rates of a cohort of randomly sampled mutants from both strains were utilized as a proxy for emerging phenotypic, and by association genetic variation. Mutant growth rates were consistently less than rates in non-mutants, an indicator of reduced fitness and the presence of mildly deleterious mutations in both the strains. In addition, the effect of these mutations on the populations as a whole varied by strain. CONCLUSIONS: Overall, this study shows that genome reductions in the MDS42 did not stabilize the chassis under metabolic stress. Over time, this could compromise the effectiveness of synthetic organisms built on chassis in environmental applications.

Single-Cell Microfluidics to Study the Effects of Genome Deletion on Bacterial Growth Behavior.

By directly monitoring single cell growth in a microfluidic platform, we interrogated genome-deletion effects in Escherichia coli strains. We compared the growth dynamics of a wild type strain with a clean genome strain, and their derived mutants at the single-cell level. A decreased average growth rate and extended average lag time were found for the clean genome strain, compared to those of the wild type strain. Direct correlation between the growth rate and lag time of individual cells showed that the clean genome population was more heterogeneous. Cell culturability (the ratio of growing cells to the sum of growing and nongrowing cells) of the clean genome population was also lower. Interestingly, after the random mutations induced by a glucose starvation treatment, for the clean genome population mutants that had survived the competition of chemostat culture, each parameter markedly improved (i.e., the average growth rate and cell culturability increased, and the lag time and heterogeneity decreased). However, this effect was not seen in the wild type strain; the wild type mutants cultured in a chemostat retained a high diversity of growth phenotypes. These results suggest that quasi-essential genes that were deleted in the clean genome might be required to retain a diversity of growth characteristics at the individual cell level under environmental stress. These observations highlight that single-cell microfluidics can reveal subtle individual cellular responses, enabling in-depth understanding of the population.

Metagenomic Sequencing Unravels Gene Fragments with Phylogenetic Signatures of O2-Tolerant NiFe Membrane-Bound Hydrogenases in Lacustrine Sediment.

Many promising hydrogen technologies utilising hydrogenase enzymes have been slowed by the fact that most hydrogenases are extremely sensitive to O2. Within the group 1 membrane-bound NiFe hydrogenase, naturally occurring tolerant enzymes do exist, and O2 tolerance has been largely attributed to changes in iron-sulphur clusters coordinated by different numbers of cysteine residues in the enzyme's small subunit. Indeed, previous work has provided a robust phylogenetic signature of O2 tolerance [1], which when combined with new sequencing technologies makes bio prospecting in nature a far more viable endeavour. However, making sense of such a vast diversity is still challenging and could be simplified if known species with O2-tolerant enzymes were annotated with information on metabolism and natural environments. Here, we utilised a bioinformatics approach to compare O2-tolerant and sensitive membrane-bound NiFe hydrogenases from 177 bacterial species with fully sequenced genomes for differences in their taxonomy, O2 requirements, and natural environment. Following this, we interrogated a metagenome from lacustrine surface sediment for novel hydrogenases via high-throughput shotgun DNA sequencing using the Illumina™ MiSeq platform. We found 44 new NiFe group 1 membrane-bound hydrogenase sequence fragments, five of which segregated with the tolerant group on the phylogenetic tree of the enzyme's small subunit, and four with the large subunit, indicating de novo O2-tolerant protein sequences that could help engineer more efficient hydrogenases.

Somatic instability of the expanded CTG triplet repeat in myotonic dystrophy type 1 is a heritable quantitative trait and modifier of disease severity.

Deciphering the contribution of genetic instability in somatic cells is critical to our understanding of many human disorders. Myotonic dystrophy type 1 (DM1) is one such disorder that is caused by the expansion of a CTG repeat that shows extremely high levels of somatic instability. This somatic instability has compromised attempts to measure intergenerational repeat dynamics and infer genotype-phenotype relationships. Using single-molecule PCR, we have characterized more than 17 000 de novo somatic mutations from a large cohort of DM1 patients. These data reveal that the estimated progenitor allele length is the major modifier of age of onset. We find no evidence for a threshold above which repeat length does not contribute toward age at onset, suggesting pathogenesis is not constrained to a simple molecular switch such as nuclear retention of the DMPK transcript or haploinsufficiency for DMPK and/or SIX5. Importantly, we also show that age at onset is further modified by the level of somatic instability; patients in whom the repeat expands more rapidly, develop the symptoms earlier. These data establish a primary role for somatic instability in DM1 severity, further highlighting it as a therapeutic target. In addition, we show that the level of instability is highly heritable, implying a role for individual-specific trans-acting genetic modifiers. Identifying these trans-acting genetic modifiers will facilitate the formulation of novel therapies that curtail the accumulation of somatic expansions and may provide clues to the role these factors play in the development of cancer, aging and inherited disease in the general population.

Association of reading disabilities with regions marked by acetylated H3 histones in KIAA0319.

Reading disabilities (RDs) have been associated with chromosome 6p with recent studies pointing to two genes, DCDC2 and KIAA0319. In this study, markers across the 6p region were tested for association with RD. Our strongest findings were for association with markers in KIAA0319, although with the opposite alleles compared with a previous study. We also found association with markers in VMP, but not with DCDC2. Current evidence indicates that differential regulation of KIAA0319 and DCDC2 contributes to RD, thus we used chromatin immunoprecipitation coupled with genomic tiling arrays (ChIP-chip) to map acetylated histones, a molecular marker for regulatory elements, across a 500 kb genomic region covering the RD locus on 6p. This approach identified several regions marked by acetylated histones that mapped near associated markers, including intron 7 of DCDC2 and the 5' region of KIAA0319. The latter is located within the 70 kb region previously associated with differential expression of KIAA0319. Interestingly, five markers associated with RD in independent studies were also located within the 2.7 kb acetylated region, and six additional associated markers, including the most significant one in this study, were located within a 22 kb haplotype block that encompassed this region. Our data indicates that this putative regulatory region is a likely site of genetic variation contributing to RD in our sample, further narrowing the candidate region.

Association of attention-deficit/hyperactivity disorder with a candidate region for reading disabilities on chromosome 6p.

BACKGROUND: Reading disabilities (RD) and attention-deficit hyperactivity/disorder (ADHD) are two common childhood disorders that co-occur by chance more often than expected. Twin studies and overlapping genetic linkage findings indicate that shared genetic factors partially contribute to this comorbidity. Linkage of ADHD to 6p, an identified RD candidate locus, has previously been reported, suggesting the possibility of a pleiotropic gene at this locus. RD has been previously associated with five genes in the region, particularly DCDC2 and KIAA0319. METHODS: To test whether these genes also contribute to ADHD, we investigated markers previously associated with RD for association with ADHD and ADHD symptoms in a sample of families with ADHD (n = 264). Markers were located in two subregions, VMP/DCDC2 and KIAA0319/TTRAP. RESULTS: Across all analyses conducted, strong evidence for association was observed in the VMP/DCDC2 region. Association was equally strong with symptoms of both inattention and hyperactivity/impulsivity, suggesting that this locus contributes to both symptom dimensions. Markers were also tested for association with measures of reading skills (word identification, decoding); however, there was virtually no overlap in the markers associated with ADHD and those associated with reading skills in this sample. CONCLUSIONS: Overall this study supports a previous linkage study of ADHD indicating a risk gene for ADHD on 6p and points to VMP or DCDC2 as the most likely candidates.

The KIAA0319-like (KIAA0319L) gene on chromosome 1p34 as a candidate for reading disabilities.

A locus on chromosome 1p34-36 (DYX8) has been linked to developmental dyslexia or reading disabilities (RD) in three independent samples. In the current study, we investigated a candidate gene KIAA0319-Like (KIAA0319L) within DYX8, as it is homologous to KIAA0319, a strong RD candidate gene on chromosome 6p (DYX2). Association was assessed by using five tagging single nucleotide polymorphisms in a sample of 291 nuclear families ascertained through a proband with reading difficulties. Evidence of association was found for a single marker (rs7523017; P=0.042) and a haplotype (P=0.031), with RD defined as a categorical trait in a subset of the sample (n=156 families) with a proband that made our criteria for RD. The same haplotype also showed evidence for association with quantitative measures of word-reading efficiency (i.e., a composite score of word identification and decoding; P=0.032) and rapid naming of objects and colors (P=0.047) when analyzed using the entire sample. Although the results from the current study are modestly significant and would not withstand a correction for multiple testing, KIAA0319L remains an intriguing positional and functional candidate for RD, especially when considered alongside the supporting evidence for its homolog KIAA0319 on chromosome 6p. Additional studies in independent samples are now required to confirm these findings.

The C2H2 zinc-finger protein SYD-9 is a putative posttranscriptional regulator for synaptic transmission.

Communication between neurons is largely achieved through chemical synapses, where neurotransmitters are released from synaptic vesicles at presynaptic terminals to activate postsynaptic cells. Exo- and endocytosis are coordinated to replenish the synaptic vesicle pool for sustained neuronal activity. We identified syd-9 (syd, synapse defective), a gene that encodes multiple C2H2 zinc-finger domain-containing proteins specifically required for synaptic function in Caenorhabditis elegans. syd-9 loss-of-function mutants exhibit locomotory defects, a diffuse distribution of synaptic proteins, and decreased synaptic transmission with unaffected neurodevelopment. syd-9 mutants share phenotypic and ultrastructural characteristics with mutants that lack synaptic proteins that are required for endocytosis. syd-9 mutants also display genetic interactions with these endocytotic mutants, suggesting that SYD-9 regulates endocytosis. SYD-9 proteins are enriched in the nuclei of both neuron and muscle cells, but their neuronal expression plays a major role in locomotion. SYD-9 isoforms display a speckle-like expression pattern that is typical of RNA-binding proteins that regulate premRNA splicing. Furthermore, syd-9 functions in parallel with unc-75 (unc, uncoordinated), the C. elegans homologue of the CELF/BrunoL family protein that regulates mRNA alternative splicing and processing, and is also required specifically for synaptic transmission. We propose that neuronal SYD-9 proteins are previously uncharacterized and specific posttranscriptional regulators of synaptic vesicle endocytosis.